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ABSTRACT 

The  objective  functions  used  in  Engineering  Optimization  are  complex  in  nature  with  many  variables  and 
constraints.  Conventional  optimization  tools  sometimes  fail  to  give  global  optima  point.  Very  popular  methods  like 
Genetic  Algorithm,  Pattern  Search,  Simulated  Annealing,  and  Gradient  Search  are  useful  methods  to  find  global  optima 
related  to  engineering  problems.  This  paper  attempts  to  use  new  non-traditional  optimization  algorithms  which  are  used  to 
find  the  minimum  cost  of  designing  welded  beam  to  obtain  global  optimum  solutions.  The  cost,  number  of  iterations  and 
the  total  elapsed  time  to  complete  the  problems  are  all  compared  using  these  ten  non-traditional  optimization  methods. 

KEYWORDS:  Welded  Beam  Design,  Pattern  Search,  Simulate  Annealing,  Pattern  Search,  GODLIKE,  Cuckoo  Search, 
Firefly  Algorithm,  Flower  Pollination,  Ant  Lion  Optimizer,  Gravitational  Search  Algorithm,  Multi- Verse  Optimizer 

INTRODUCTION 

Welding  can  be  defined  as  a process  of  joining  metallic  parts  by  heating  to  a suitable  temperature  with  or  without 
the  application  of  pressure.  Welding  is  an  economical  and  efficient  method  for  obtaining  a permanent  joint  of  metallic 
parts.  There  are  two  distinct  applications  of  welded  joints. 

• A welded  joint  can  be  used  as  a substitute  for  a riveted  joint 

• A welded  structure  can  be  used  as  an  alternative  method  for  casting  or  forging. 

Welding  process  are  broadly  classified  with  the  following  two  groups 

• Welding  process  that  use  heat  alone  to  join  the  two  parts 

• Welding  process  that  use  a combination  of  heat  and  pressure  to  join  the  two  parts.  (Bhandari.  V.B) 

The  welding  process  that  uses  heat  alone  is  called  fusion  welding  process.  In  this  method  the  parts  to  be  joined  are 
held  in  position  and  molten  metal  is  supplied  to  the  joint.  The  molten  metal  can  come  either  from  the  parts  themselves 
called  Parent  metal  or  external  filler  metal  is  supplied  to  the  joint.  The  joining  surface  of  two  parts  becomes  plastic  or  even 
molten  under  the  action  of  heat.  When  the  joint  solidifies,  the  two  parts  fuse  into  a single  unit. 

A beam  is  a member  subjected  to  loads  applied  transverse  to  the  long  dimension,  causing  the  member  to  bend. 
Beams  are  frequently  classified  on  the  basis  of  supports  or  reactions.  A beam  supported  by  pins,  rollers  or  smooth  surfaces 
at  the  ends  is  called  simple  beam.  A beam  support  will  develop  a reaction  normal  to  the  beam  but  will  not  produce  a 
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moment  at  the  reactions.  If  either  or  both  ends  of  a beam  projects  beyond  the  supporters,  it  is  called  simple  beam  with 
overhang. 

Description 

A beam  A is  to  be  welded  to  a rigid  member  B.  The  welded  beam  is  to  consist  of  1010  steel  and  is  to  support  a 
force  P of  6000  lb.  The  dimensions  of  the  beam  are  to  be  selected  so  that  the  system  cost  is  minimized.  A schematic  of  the 
system  is  shown  in  the  figure  1 . 


Figure  1:  Welded  Beam 

The  diagram  illustrates  a rigid  member  (A)  welded  onto  a beam.  A load  is  applied  to  the  end  of  the  member.  The 
beam  is  to  be  optimized  for  minimum  cost  by  varying  the  weld  and  member  dimensions 
X = (. h , l,  t,  b)  = (Xj , x2 , x3 , x4  ) 

This  includes  limits  of  the  shear  stress,  bending  stress,  buckling  load  and  end  deflection.  The  variables  and  x2 
are  integer  multiples  of  0.0625  inch,  but  for  this  application  are  assumed  continuous.  (Amie  Mesari,  2012)  (Hong-Shaung 
Li,  Siu-Hui  Au,  2010) 

PARAMETERS: 

(Welded  Beam  Design  Optimization  ) 

Young’s  Modulus  (psi) 

E =30xl06  psi 

Shearing  modulus  for  the  beam  material  (psi) 

G = 12xl06  psi 

Overhang  length  of  the  member  (inch) 

L =14  in 

Design  stress  of  the  weld  (psi) 

T = 13600  psi 

max 

Design  normal  stress  for  the  beam  material  (psi) 

o =30000  psi 

max 


Impact  Factor  (JCC):  2.0346 


NAAS  Rating:  3.19 
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Maximum  deflection  (inch) 

8 =0.25  in 

max 

Load  (lb) 

P =6000  lb 

COST  FUNCTION: 

(Welded  Beam  Design  Optimization  ) 

The  performance  index  appropriate  to  this  design  is  the  cost  of  a weld  assembly.  The  major  cost  components  of 
such  an  assembly  are  (1)  set  up  labour  cost,  (2)  welding  labour  cost,  and  (3)  material  cost: 

f(X ) =CQ  +Cj  +c2 


where  f(X)  = cost  function 


C, 


C„ 


= set  up  cost 
= welding  labor  cost 
= material  cost 


Setup  Cost  r 

L0 


The  company  has  chosen  to  make  this  component  a weldment,  because  of  the  existence  of  a welding  assembly 
line.  Furthermore,  assume  that  fixtures  for  setup  and  holding  of  the  bar  during  welding  are  readily  available.  The  cost  q 

can  therefore  be  ignored  in  this  particular  total-cost  model. 

Welding  Labour  Cost  c ^ 


Assume  that  the  welding  will  be  done  by  machine  at  a total  cost  of  $10/hr  (including  operating  and  maintenance 
expense).  Furthermore,  suppose  that  the  machine  can  lay  down  a cubic  inch  of  weld  in  6 min.  The  labour  cost  is  then 


C,  =|  10 


hr 


1 hr 
60  min 


f 


\ 


min 

v67T 

V in 


V =1 

w 


Where  V = weld  volume,  in3. 

w 


Material  Cost , 


CL  =C,  V +C.VD 
2 3 w 4 B 


where 


= cost  per  volume  per  weld  material,  $/in  = (0.37)(0.283) 
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£ = cost  per  volume  of  bar  stock,  $/in  ’ = (0.17)(0.283) 


V = volume  of  bar  A,  in  . 
B 


From  the  geometry. 


V = 2|  -h2l)  = h2  l 

w 1 2 J 


Volume  of  the  weld  material  (inch3) 


Vweld  X1 X2 


V =tb(L  + l)  , 


Volume  of  bar  (inch3) 


Vbar=Xf*,L  + X2> 


Therefore  the  cost  function  become 


f(X ) =h2I  +C3h2  l + C 4tb(L  + l) 


Or,  in  terms  of  the  x variables. 


f(X)  = (l  + C3  )x2  x2  + C4  x3x4  (14.0  + x2 ) 


f(X)  =1.10471  X|  x2  + 0.048 1 1 x^  (14.0  + x2 ) 

EN GINEERIN GRELATIONSHIP : (Hasancebi.O,  2012) 

To  complete  the  model,  it  is  necessary  to  define  the  important  stress  states. 
Weld  stress 


t(X)  = 3Jt2  +2t1t2  cos  e + T2  Where  COS<9=  — 

Weld  stress  has  two  components  T|  and  T 2 , where  T ^ the  primary  stress  acting  over  the  weld  throat  area  and  X 2 
is  a secondary  torsional  stress. 


Primary  stress  acting  over  the  weld  throat 


1 JTx^x2 


Impact  Factor  (JCC):  2.0346 


NAAS  Rating:  3.19 
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Secondary  torsion  stress 

MR 

r2-T 

Moment  of  P about  centre  of  gravity  of  weld  setup 

= 84000  + 3000 


M =P 


R =. 


r V 

L + — 

2 

V V 


- + 


V*3 


Polar  moment  of  inertia  of  weld 


J = 2\ 


V2  XyX^ 


12 


- + 


"vV2 


Bar  bending  stress:  The  maximum  bending  stress  can  be  shown  to  be  equal  to 

6 PL  = 504000 


<r{X)  = 


2 2 
X4X3  X4X3 


Bar  deflection.  To  calculate  the  deflection,  assume  the  bar  to  be  a cantilever  of  length 

4 pi}  _ 2.1952 


S(x) 


E x ^ x ^ 


The  buckling  load  is  approximated  by 


2 6 

X3XA 

4.013£iM-JL  ( 3 

36 


P(X)  = 


l-: 


2 L V 4G 


= 64746.022(1  -0.028234 jc  )jc  jcJ 


CONSTRAINT: 

The  shear  stress  at  the  beam  support  location  cannot  exceed  the  maximum  allowable  for  the  material 

r(X)  < r 

max 
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T(X)  < 13600 

The  normal  bending  stress  at  the  beam  support  location  cannot  exceed  the  maximum  yield  strength  for  the 


material 


ct(X)  < cr 


max 


504000 
x x2 
4 3 


< 30000 


16.8  < x^x^ 


The  weld  thickness  is  less  than  the  material  thickness 


xl  ^ x4 


+c^x^x^{L  + x^)  < 5.0 


0. 1047 Lc“jc2  + 0.048 lh3*4 (14.0  + x2)  < 5.0 


The  defined  minimum  must  be  less  than  the  weld  thickness 

0.125  < x{ 

The  deflection  cannot  exceed  the  maximum  deflection 


S(X)  < S 


max 


2.1952 

*3*4 


< 0.25 


1.7808  < x^x4 


The  applied  load  is  less  than  the  buckling  load 


P < Pc(X) 


6000  < 64746.022(1 -0.028234x3)x3x^ 


0.0926697859  < (1  -0.028234x3)x3x^ 


Size  constraints 


Impact  Factor  (JCC):  2.0346 


NAAS  Rating:  3.19 
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x i >0.1 

x2>0.1 

*3  -10 
x4  <2 

Hence  the  minimum  cost  and  optimized  dimensions  are 
Solution  = optimization  minimize  (f(X),  cons,  bounds) 

MATHEMATICAL  FORMULATION  (Janga  Reddy.M,  D.Nagesh  Kumar,  2007)_(Ali  Riza  Yildiz,  2008) 

The  mathematical  formulation  of  the  objective  function  f(X)  which  is  the  total  fabrication  cost  mainly 
comprised  of  the  set-up,  welding  labour  and  material  cost  is  as  follows 

The  objective  is  to  minimize  the  cost  of  the  welded  beam  design  problem. 

Minimize  f(X)  = 1 .10471  +0.0481 1 x^x^ (14.0  + ) 


Subject  to 


T(X)  < T 


max 


<x(X)  < <T 


max 


*1  - x4 


0.10471x^X2 +0.0481  lx3x4 (14.0  + x2)  < 5.0 


0.125  < x 


1 


S(X)  < s 


P < P (X) 
c 


Where 


r(X)  = 


f „ X 

t}  + 2 T, 

x2 

1 1 2 

2 R 

v J 

+ Tt. 


T,  = 


1 a/2~ \ 
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MR 


Size  constraints 


0.1<x1<2 
0.1  <x2<  10 
0.1  < x3  < 10 
0.1  <x4<  2 


Impact  Factor  (JCC):  2.0346 


NAAS  Rating:  3.19 
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COMPARATIVE  RESULTS 


Trial'Noa 

PSOa 

SAo 

PSo 

GLo 

Cuckooa 

FFa 

FPa 

ALOa 

GSAa 

MYOa 

rhicknessoftherveldXi: 

0.206412a 

0.165306; 

0.55; 

0.204164; 

0.20573; 

0.214698; 

0.205729; 

0.177859; 

0.219556: 

0.199033: 

Lengthoftheweldx:: 

3.528353; 

5.294754; 

2.028512; 

3.565391; 

3.519497; 

3.655292; 

3.519502; 

4.393466: 

4.728342: 

3.652944: 

Lengthofthematerialx]: 

8.988437; 

8.872164; 

4.798717; 

9.05924; 

9.036624; 

8.507188; 

9.036628: 

9.065462: 

8.50097; 

9.114448: 

Thickness  ofthe  material  Xf: 

0.208052; 

0.217625; 

0.729557; 

0.206216; 

0.20573; 

0.234477; 

0.20573: 

0.20559: 

0.271548: 

0.205478: 

"ostn 

1.742326o 

1.939196; 

3.37756; 

1.7428; 

1.731527; 

1.864164; 

1.731528: 

1.796793: 

2.295076: 

1.749834: 

iteration; 

178; 

4666; 

5a 

66573; 

100000; 

20000; 

2000; 

1000; 

1000: 

1000: 

[rime; 

1.235504a 

4.38926; 

0.746971; 

8.338414; 

12.45931; 

8.083191; 

7.041082: 

59.04323: 

23.97281: 

8.962193: 
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Impact  Factor  (JCC):  2.0346 


NAAS  Rating:  3.19 


Cost  Minimization  of  Welded  Beam  Design  Problem  Using 
PSO,  SA,  PS,  Godlike,  Cuckoo,  FF,  FP,  ALO,  GSA  and  MVO 


11 


V) 

C 

o 


ru 

i- 

<D 


120000 

100000 

80000 

60000 

40000 

20000 

0 


t 

; 

_L 

o < 

IS)  </i 
Q- 


ID 


o 

o 

u 

3 

u 


o < o 

j w s 
< ID  ^ 


Methods 


<u 

£ 

4-» 

■a 

<u 

wi 

Cl 

_ro 

<u 


70 

60 

50 

40 

30 

20 

10 

0 


ill 


£ 


o < 

(/l  l/l 
Q- 


(/I 

O- 


ID 


O u_ 

cl  o 

< 

O LJ- 

u-  Zj 

LD 

U 

< 

CD 

u 

O 

> 


methods 


RESULTS  AND  DISCUSSION 


With  the  two  extreme  values  of  the  parameters  the  optimization  is  carried  out  with  different  solvers.  As  they  are 
stochastic  type  the  results  may  vary  from  trial  to  trial.  So  the  problem  is  made  to  run  for  20  trials.  (Elbeltagi.E.,  Tarek 
Hegazy.I.,  Grierson  D.,  2005)  And  an  average  of  all  trials  are  taken  as  a final  value  of  the  parameter  by  the  solver.  The 
solvers  are  compared  with  three  different  criteria. 

1.  Consistency: 

The  cost  is  consistent  in  pattern  search  (3.37756) 

2.  Minimum  run  time: 

For  minimum  run  time  of  the  problem  we  have  PS  (0.746971  seconds),  PSO  (1.235504  seconds). 

3.  Minimum  Evaluation: 

This  Criterion  will  determine  the  effectiveness  of  the  algorithm.  From  the  table  we  see  that  the  PS  and  PSO 
algorithm  have  minimum  evaluation  of  5 and  178  respectively. 

4.  The  Simplicity  of  Algorithm: 

Of  all  the  algorithms.  Pattern  Search  algorithm  is  the  most  simplest  followed  by  Particle  Swarm  Optimization. 
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Thus  it  is  seen  that  the  PS  solver  satisfies  all  the  criteria.  Even  though  the  pattern  search  satisfies  all  the  above 
criteria,  the  cost  becomes  maximum  whereas  the  cost  in  PSO  is  1.7423.  Therefore  the  particle  swarm  optimization  has  the 
minimum  cost  with  time  1.23  seconds  and  178  iteration.  So  the  appropriate  algorithm  for  welded  beam  design  is  suggested 
as  Particle  Swarm  Optimization. 

CONCLUSIONS 

In  the  present  study  the  PSO  algorithm  is  proposed  as  a simple  and  efficient  optimization  technique  for  handling 
welded  beam  design  problem.  PSO  algorithm  is  a population  based  technique  which  follows  a stochastic  iterative 
procedure  to  locate  the  optimum  or  a reasonably  near-  optimum  solution  for  the  welded  beam  design  optimization. 
Performance  evaluation  of  the  PSO  algorithm  through  welded  beam  design  optimization  reveals  the  efficiency  of  this 
technique  in  solving  practical  optimization  problems.  Although  in  the  present  study  the  PSO  algorithm  is  utilized  only  for 
solving  welded  beam  design  optimization  problem,  it  can  be  easily  employed  for  solving  other  types  of  optimization 
problems  as  well. 
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APPENDIES 

Tables  having  the  option  set  and  stopping  criteria  for  the  optimization  methods 


Methods 

PSO 

SA 

PS 

GL 

CUCKOO 

FF 

FP 

ALO 

GSA 

MVO 

Option  it 

t 

MasJSmsuwa 

=:oo 

Limit3  cc 

Average  change 

in  fitness 

value3 10-5 

Function 

Tolerance:  1(H 

Cognitive 

Attraction=0.5 

Population 

Size=40 

Social 

Attraction3!  .25 

Initial 

Temperature:  100 
Annealing 
Function:  Fast 
Annealing 

taBwatew 
Interval:  100 
Time  limit  :cc 

Evaluation's* 

SSUfifvariables. 

Function 
Tolerance:  10-5 
Obj  ective 
Limit:10-5 

Poll 

positive  basis 
2N 

Initial  Mesh 

Size:l 

Expansion 

Function:  2 

Contraction 

Factor:0.5 

Mesh 

Tolerence:  10-5 
Evaluation's 

STJisufif 

variables. 

Max. Iteration:! 

oo*3swiX 

variables. 

MasJsae. 

Limit:co 
Function 
Tolerance:  10-5 

Ma&fwi 

Evaluations 

10--' 

iUsUKMB 

snr-° 

Jfi&Jssati 

star- 

Total 

Iterations= 

15 

Functions 

Tolerance= 

10- 

Jtto^aa 

Evaluations 

=10-5 

MasJisaaa 

as=:o 

Functions 

Tolerance= 

io-« 

Limit3  oc 

Jttofaa 

Evaluations 

=10-5 

MasJisaaa 

OS=- o 
Functions 
Tolerance= 
10-« 

MasJsat 

Limit3  oo 

MasJjw 

Evaluations 

=10-5 

as=:o 

Functions 

Tolerance3 

10-5 

MaxC&ae. 

Limit3  cc 

M&sbfaa 

Evaluations=l 
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